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Application  of  single-beam  reflective  laser  optical  interferometry  for  oil  films  and  droplets  in  water  de¬ 
tection  and  characterization  is  discussed  Oil  films  can  be  detected  by  the  appearance  of  characteristic 
interference  patterns  Analytical  expressions  describing  intensity  distribution  in  these  interference  pat¬ 
terns  allow  determination  of  oil  film  thickness,  size  of  oil  droplets,  and  distance  to  the  oil  film  from  the 
observation  plane.  Results  from  these  analyses  indicate  that  oil  spill  aging  and  breakup  can  be  monitored 
in  real  time  by  analyzing  time-dependent  holographic  fringe  patterns.  Interferometric  methods  of  oil  spill 
detection  and  characterization  can  be  automated  using  digital  holography  with  three-dimensional  recon¬ 
struction  of  the  time-changing  oil  spill  topography.  In  this  effort,  the  interferometric  methods  were  ap¬ 
plied  to  samples  from  Chevron  oil  and  British  Petroleum  MC252  oil  obtained  during  the  Deep  Water 
Horizon  oil  spill  in  the  Gulf  of  Mexico.  ©  2011  Optical  Society  of  America 
OCIS  codes:  090.0090,  090  2880,  090.5694. 


1.  Introduction 

Optical  interferometry  is  a  practical  and  popular 
method  of  detection  and  characterization  of  oil  films 
and  emulsified  liquids.  In  this  paper,  we  consider  re¬ 
flection  geometry  and  briefly  discuss  different  meth¬ 
ods  of  optical  interferometry  with  a  focus  on  the 
detection  and  quantification  of  the  thickness  of  oil 
spills  on  the  water.  Theoretical  solution  of  the  Max¬ 
well  wave  equations  [1]  for  reflection  of  electromag¬ 
netic  waves  from  the  oil  films  and  droplets  requires 
proper  boundary  and  initial  conditions.  Backscatter- 
ing  from  emulsified  liquids  was  investigated  in  [2], 
and  coherent  backscattering  measurement  was  ap¬ 
plied  for  the  study  of  microstructures  in  emulsified 
liquids.  An  enhanced  component  toward  the  backscat¬ 
tering  direction  emerges  superimposed  on  the  diffuse 
light  scattering  from  a  sample  with  randomly 
dispersed  scatterers,  giving  evidence  of  the  weak  lo¬ 
calization  of  light  in  random  systems.  Accurate  mea- 
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surement  of  the  coherent  backscattering  peak  gives 
us  information  on  the  optical  transport  mean  free 
path  in  the  sample,  which  leads  to  an  estimation  of 
the  droplet  size  in  emulsion  A  new  measurement  sys¬ 
tem  was  developed  and  clear  backscattering  peaks 
were  successfully  observed  for  various  kinds  of  emul¬ 
sions.  The  key  technique  of  the  current  experiment 
was  to  decrease  the  laser  speckle  by  stirring  the  li¬ 
quid.  The  result  obtained  for  heptane/AOT  [sodium 
bis(2-ethylhexyl)]/water  emulsion  was  qualitatively 
related  to  the  droplet  size,  which  could  be  controlled 
by  changing  the  amount  of  added  salt.  The  effect  of 
excluded  volume  on  the  scattering  mean  free  path 
was  discussed. 

The  reflectance  of  oil  emulsion  grows  when  wave¬ 
length  increases  with  decreasing  size  of  oil  droplets 
in  emulsion  [3].  Another  possibility  for  oil-film  thick¬ 
ness  determination  is  white-light  interferometry 
[4,5],  but  it  needs  a  more  complex  algorithm  for  de¬ 
ciphering  the  interference  patterns  in  comparison 
with  laser  monochromatic  interferometry. 

The  laser  photophoretic  behavior  of  an  oil  droplet 
in  water  emulsions  was  examined  by  radiating  a 
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single-mode  Nd:YAG  laser  (1064nm)  [6].  In  the  [6] 
emulsion  systems  examined,  the  photophoretic  ve¬ 
locity  (/vm  s'1)  of  an  organic  droplet  was  found  to  be 
linearly  proportional  to  both  the  laser  power  and  the 
radius  of  the  droplet.  Furthermore,  for  droplets  of  the 
same  size,  the  photophoretic  mobility  (/^ms“1W"1), 
which  means  the  photophoretic  velocity  normalized 
with  the  laser  power,  increased  as  the  refractive  in¬ 
dex  of  the  droplet  increased.  The  dependence  of  the 
photophoretic  mobility  was  discussed  on  the  basis  of 
a  simple  radiation  theory  in  which  the  scattering 
force  of  a  laser  beam  on  a  spherical  droplet  was 
formulated. 

In  this  paper,  we  utilize  samples  from  Chevron  oil 
and  British  Petroleum  (BP)  light  crude  MC252  oil. 
We  consider  laser  single-wavelength  and  single¬ 
beam  interferometry  as  the  most  simple  and  reliable 
method  of  oil-film  slick  detection  and  quantification 
of  the  oil-film  thickness  and  oil  droplets.  Reflected 
beams  from  the  oil  films  and  droplets  on  the  water 
are  described  by  a  holographic  approach  as  the  inter¬ 
ference  of  the  signal  wave  reflected  from  the  top  of 
the  oil  spot  with  the  reference  beam  reflected  from 
the  oil-water  interface  (for  films)  or  from  the  water 
surface  (for  droplets  and  islands). 

2.  Experimental  Setup  and  Modeling 


(a)  <b) 

Fig.  2.  (Color  online)  Typical  interference  patterns  on  reflection 
from  the  oil-on-water  film  under  near-normal  illumination  by  laser. 

mation  from  the  interference  pattern,  proper  model¬ 
ing  of  fringe  pattern  formation  is  needed. 

B.  Modeling  of  Reflection  from  the  Oil-on-Water  Film 

To  find  the  reflection  coefficient  of  the  light  wave 
with  amplitude  Eo  from  the  oil-on-water  film,  we  will 
consider  a  three-layer  system  (Fig.  3)  air — 1,  oil  film 
— 2,  and  water — 3,  with  refractive  indices  n2,  and 
n3.  The  total  reflected  wave  amplitude  Er  will  be  re¬ 
presented  as  a  sum  of  the  reflected  wave  from  the 
front  surface  Ef  and  from  the  back  surface  Eb : 

Er  =  Ef  +  Eb  =  (r  \2E$  +  (1) 


A.  Experimental  Results 

A  simple  scheme  of  holographic  interferometry  is 
shown  in  Fig.  1.  A  beam  from  an  He-Ne  laser  (wave¬ 
length  0.632 /im),  expanded  by  a  lens  (F  =  5cm), 
illuminates  a  sample  of  oil-on-water  film  at  the 
near-vertical  direction.  The  diameter  of  illuminated 
area  at  the  sample  is  about  3  cm.  Optical  beams  re¬ 
flected  from  the  top  and  bottom  surfaces  of  the  oil 
film  form  interference  patterns  on  the  screen  (which 
was  placed  at  a  6  m  distance  from  the  oil  film).  This 
interference  pattern  may  be  captured  by  the  CCD 
camera  and  digitally  processed.  Examples  of  typical 
interferogram  are  shown  in  Fig.  2;  photos  were  taken 
from  the  white  diffusion  screen.  Because  of  differ¬ 
ences  in  reflectance,  thick  oil  film  is  seen  as  bright 
spots,  while  water  appears  as  dark  spots  in  the  pic¬ 
ture.  The  geometry  of  fringes  reflects  the  thickness 
gradient  of  oil  film  and  contains  information  about 
the  oil-film  topography.  To  recover  topography  infor- 

I  ascr  Mirror 


Fig.  1.  Scheme  of  oil-on-water  film  interferometry;  distance 
between  oil  film  and  screen  ~6m. 


Here  r12  and  r23  are  reflection  coefficients  from  air- 
oil  and  oil-water,  respectively,  tl2  and  f21  are  trans¬ 
mission  coefficients  from  air-oil  and  from  oil-air,  k  is 
the  wavenumber  k  =  2;r/A,  /  is  the  distance  along  the 
propagation  direction,  and  O  is  the  phase  shift 
between  these  two  wave  due  to  the  optical  path 
difference: 

O  =  2 kh\nl  -  (nx  sinQ)2)1/2.  (2) 

In  derivation  of  Eq.  (2),  we  ignore  multiple  reflec¬ 
tions,  which  is  justified  for  the  small  reflection  coef¬ 
ficients  of  the  oil-water  interface.  For  the  slanted 
geometry  of  observations  (not  strictly  vertical  illu¬ 
mination),  amplitude  transmission/reflection  coeffi¬ 
cients  depend  on  polarization  that  may  be  used  for 
improvement  of  detection  sensitivity. 

For  s-polarization  [horizontal,  perpendicular  to  the 
plane  of  incidence,  or  transverse  electric  (TE)],  the 
coefficients  for  the  air-water  interface  are  [1] 


Fig.  3.  Scheme  of  reflection  of  the  wave  with  amplitude  Eo  from 
the  oil  film  of  thickness  h  w  ith  refractive  index  n2  (refractive  in¬ 
dices  of  air  and  water  are  n,  and  n3,  respectively). 
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rii  cos  8  -  n2  cos  8, 
ri2  nx  cos  8  +  n2  cos  8, 
2  n*  c 


(3) 


Here  8,  is  the  transmission  angle,  which  is  related  to 
the  incidence  angle  by 


cos  8,  = 


sin8  j2 


It  follows  from  Eq.  (3)  that,  for  n2  >  n\  (as  in  our  case 
of  the  air-water  interface),  the  reflection  coefficient 
is  negative,  which  is  the  equivalent  of  a  ;r-phase  shift 
during  reflection  of  s-polarized  light. 

For  p-polarization  [vertical,  parallel  to  the  plane  of 
incidence,  transverse  magnetic  (TM)],  the  coeffi¬ 
cients  will  take  the  form 

n\  cos  8r  -  n2cos8 

r12=  — - K- - - - k  •  4) 

nx  cos  8,  +  rc2cos^ 

For  the  transmission  coefficient,  we  can  get  the 
expression 


reference  wave  created  by  reflection  of  the  plane 
wave  from  the  water  surface.  The  signal  wave  will 
be  produced  by  reflection  of  part  of  the  laser  wave 
from  the  curved  surface  of  the  oil  droplet.  Interfer¬ 
ence  of  these  copropagating  waves  produces  inter¬ 
ference  patterns  that  depend  on  the  oil  droplet 
thickness. 

Scattering  of  light  by  spherical  particles  has  been 
intensively  studied  by  Mie  theory  [10,11]  and  later 
was  extended  to  nonspherical  particles  [11,12].  The 
special  case  of  the  forward  scattering  in  the  far  zone 
may  be  also  described  as  single-beam  (or  in-line) 
holography. 

From  Mie  theory  [10],  the  backward  scattered 
electric  field  is  given  by 

E,=e-ikr-h*SE0/(ikr).  (7) 

Here,  k  is  the  wave  number,  S  is  the  scattering  am¬ 
plitude,  and  E0  is  the  incident  plain  wave.  In  the 
point  x ,  y,  z  in  the  backward  direction,  since  x,  y 
z  we  have 


Z  4 


x2  +y2 
2 2  ' 


(8) 


.  cos  8 

*12  =  (l  +  r12) - K“*  (5) 

COS  8, 

For  this  TM  polarization,  there  exists  a  Brewster  an¬ 
gle  of  no  reflection,  defined  by  equating  the  reflection 
coefficient  in  Eq.  (3)  to  zero. 

The  general  approach  for  analysis  of  interference 
patterns  is  to  find  the  first  phase  shift  O  using  the 
expressions  [7,8] 


<t>  =  tan 

mux))  =  f-> 


where 

f(A)f(/w>  . 


(6) 


Superposition  of  the  incident  (reference)  and 
scattered  (signal)  waves  give  total  amplitude 

E  =  Es  +  E0  =  E0[  1  4  e )/2zS / (ikz)\ .  (9) 

For  intensity  I  in  approximation  of  Es  «:  E0}  we  will 
get 

Re[-jSe-‘*(x24v2)/2z]j,  (10) 

where  Re  =  real  part 

For  radius  of  the  scattering  particle  a  » 
(laser  wavelength),  the  scattering  amplitude  is,  for 
perpendicular  polarization  [13,14], 


and  H(I(x))  is  the  Hilbert  transform  of  the  intensity 
signal,  F  represents  the  Fourier  transform,  and  F  1 
is  the  inverse  Fourier  transform.  Once  the  phase  has 
been  determined,  the  oil-film  height  can  be  found 
from  Eq.  (2).  As  it  was  mentioned  in  [9],  comparison 
of  TE  and  TM  polarizations  for  the  film  interferome¬ 
try  show  that  TM  polarization  gives  better  visibility 
and  so  is  more  sensitive  for  oil-film  characterization. 

C.  Reflection  from  the  Oil  Droplets 

For  those  types  of  oil  whose  refractive  indices  are  not 
very  different  from  that  of  water,  thin-film  interfer¬ 
ence  will  give  low-contrast  patterns  that  are  difficult 
to  analyze.  In  such  cases,  interference  patterns  aris¬ 
ing  from  the  reflection  from  oil  droplets  and  islands 
will  be  more  pronounced.  Modeling  of  interference 
patterns  in  this  case  can  be  done  in  a  way  that  is  si¬ 
milar  to  wavefront  division  interferometry  with  a 


S(8)  =  /;2  [Ji(y)/r  + 


ip'Z(fi)  J0(r) 


,  sin  6 

£(<$)  - - cos  S  +  i 


d2  () 

cos  6 


-  sin  S 


»=\fr2+p2, 

Y  —  2ft  sin  8/2, 
p'  =  2(n'-l)fty 


n  2n 

p  T«. 

P  =  2(/i  -  1  )ft, 


(11) 


where  n  is  the  refractive  index  of  spherical  particles 
(oil  droplet)  n  =  n'-in"f  Jx  and  J0  are  the  Bessel 
functions  of  the  first  and  zero  orders,  0  is  the  polar 
angle  and  a  is  the  radius  of  curvature  for  the  sphe¬ 
rical  particle. 

The  first  term  i n  Eq.  ( 1 1 )  is  the  d iffracted  term .  The 
second  term  includes  the  refracted  light  term  (in 
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geometric  optics  approximation)  and  depends  on  the 
refractive  index.  The  last  term  is  more  pronounced  for 
smaller  particles  and  is  due  to  interaction  between 
the  refracted  and  diffracted  components. 

From  the  analysis  of  the  extrema  of  the  interfer¬ 
ence  pattern  [10],  the  radius  of  the  particle  a  and 
the  distance  of  the  particle  to  the  screen  z  can  be 
found. 

For  illustration,  we  will  consider  analysis  of  the 
interference  pattern  when  only  the  first,  diffraction 
term  is  important. 

To  find  the  size  of  the  diffracting  droplet,  we  need 
to  equal  the  Bessel  argument  kar\/z  to  the  values 
that  reduce  it  to  zero  (here  rx  =  (x2  4-y2)1/2). 

As  an  example  of  both  radial  and  density  plots, 
extrema  (Fig.  4)  of  sin  (at  ~0.16cm)  and  Bessel  func¬ 
tions  (at  ~  0.6  cm)  are  easily  noticeable,  which  allows 
us  to  find  distance  z  =  800  cm  and  droplet  size 
a  =  0.5  mm. 

3.  Kinetics  of  Oii-Fiim  Aging  and  Dispersion  by 
Corexlt  9500 

Application  of  dispersant  Corexit  9500  to  the  oil  sam¬ 
ples  allows  monitoring  of  the  oil  breakup  and  aging 
[15].  Figure  5  depicts  a  sequence  of  photos  that  shows 
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(a) 


Fig.  4.  (Color  online)  (a)  Density  plot  and  (b)  radial  dependence  of 
the  hologram  intensity  for  parameters  a  —  0  5  mm,  z  —  800cm, 
wavelength  0  633/^m). 


(a)  (b) 


Fig.  5.  (Color  online)  Kinetics  of  oil-film  aging  after  applying  dis¬ 
persant  Corexit  9500  (top  left  to  bottom  right);  photos  were  taken 
at  3  s  intervals  (diameter  of  the  illuminated  oil-film  area  was  about 
3  cm).  The  last  (bottom)  picture  shows  breakup  of  a  big  oil  island 
with  the  formation  of  droplets. 

the  kinetics  of  the  oil  spot  thinning  and  breaking  up 
upon  application  of  the  liquid  dispersant  to  the  film. 
All  photos  were  taken  from  the  images  projected  on 
the  white  diffusion-scattering  screen  Figure  5  shows 
that  only  the  relatively  smooth  part  of  the  kinetic, 
initial  stage  of  oil  dispersion  was  much  faster  (about 
1  s)  and  was  followed  by  the  abrupt  spatial  separa¬ 
tion  of  fragments.  After  this  (fast)  initial  stage, 
smoother  changes  with  relatively  small  separation 
were  captured  by  the  interference  method,  as  shown 
in  Fig.  5. 

For  observation  of  the  initial  (fast)  stage  of  the  oil- 
film  dispersion,  modification  of  the  interferometric 
experiment  is  needed,  including  illumination  with 
a  larger-diameter  laser  spot  to  allow  capture  of  the 
diverging  oil-film  fragments. 

4.  Discussion  and  Conclusion 

We  have  discussed  application  of  single-beam  laser 
optical  interferometry  for  detection  and  characteri¬ 
zation  of  oil  films,  islands,  and  droplets  on  water. 
Detection  of  oil  films  can  be  done  by  observation  of 
characteristic  interference  patterns.  Analytical  ex¬ 
pressions  describing  intensity  distribution  in  these 
interference  patterns  allows  finding  the  thickness 
of  the  oil  films,  the  size  of  oil  droplets,  and  the  dis¬ 
tance  to  the  oil  film  from  the  observation  plane. 

Estimates  of  the  oil-film  thickness  can  be  accom¬ 
plished  by  counting  fringes  on  the  interferogram 
[Fig.  2(b)]  at  the  edge  of  the  oil  spot.  By  multiplying 
the  number  of  fringes  (about  six)  on  the  laser  wave¬ 
length  (0.633/un),  we  get  a  film  thickness  value  of 
3.8  /^m.  By  modeling  of  the  interferogram  of  a  droplet, 
we  show  that  both  droplet  size  and  its  distance  to  the 
observation  plane  can  be  found  from  the  analysis  of 
the  interferogram  fringes.  In  addition,  it  is  possible 
to  observe  the  kinetics  of  oil  slick  aging  by  monitor¬ 
ing  the  time  dependence  of  fringes.  Interferometric 
methods  of  oil  spill  detection  and  characterization 
can  be  automated  by  using  a  decision-making  digi¬ 
tal  approach  with  proper  programming  [9]  that  in¬ 
cludes  three-dimensional  reconstruction  of  the  oil 
topography. 
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